Numerical Study of the Spin-Flop Transition in Anisotropic Spin-1/2 

Antiferromagnets 
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Magnetization processes of the spin-1/2 antiferromagnetic XXZ model in two and three spatial 
dimensions are studied using quantum Monte Carlo method based on stochastic series expansions. 
Recently developed operator-loop algorithm enables us to show a clear evidence of the first-order 
phase transition in the presence of an external magnetic field. Phase diagrams of closely related sys- 
tems, hard core bosons with nearest-neighbor repulsions, are also discussed focusing on possibilities 
of phase-separated and supersolid phases. 
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There have been a lot of interests on magnetic proper- 
ties of anisotropic quantum antiferromagnets since Neel 
first predicted the first-order phase transition in the pres- 
ence of an external magnetic field [Q. One of the sim- 
plest models for anisotropic antiferromagnets in an exter- 
nal magnetic field is described by the following spin-1/2 
XXZ model Hamiltonian: 

H^jJ2 {SiSf + SfSf + A5f 5/) ~hJ2Si, (1) 

(iJ> i 

where S" is the a(= x,y,z) component of the spin- 
1/2 operator at site i, h is an external magnetic field 
applied in the z-direction, and runs over all the 

nearest-neighbor pairs of spins at sites i and j. J(> 0) 
is antiferromagnetic coupling and A(> 0) is anisotropic 
constant. Mean-field calculations of the spin-1/2 XXZ 
model supporting Neel's prediction, found the first- 
order phase transition from the Neel ordered state to the 
spin-fiopping state with increasing the magnetic field h. 
On the contrary to these studies, it is known from Bcthc 
ansatz solution that the one-dimensional (ID) spin-1/2 
XXZ model shows a second-order transition in the pres- 
ence of the external magnetic field Q . The discrepancy 
is assigned to the inadequacy of treating quantum fluctu- 
ations by the mean-field theories. It is thus important to 
use an unbiased numerical method for understanding the 
correct nature of the magnetization process even for the 
simplest systems such as one given by Eq. (|^) in higher 
spatial dimensions since there exist no analitically exact 
solutions. This is preciely one of our purposes for this 
paper [||. 

Another importance of studing the spin-1/2 XXZ 
model defined by Eq. (|l|) comes from the fact that the 
model is mapped onto a system of hard core bosons with 
nearest-neighbor repulsions described by the Hamilto- 



(iJ) 



V ^ ninj 
(iJ> 



(2) 



with t = J/2, V = J A, and fi = h + zJA/2 (2: coor- 
dination number) Q. Here Cj is a creation operator of 



hard core boson at site i and Ui = CjCi. The total mag- 
netization Mz = J2i^i t^us relates to J2i ("i ~ 1/2) in 
the boson model. The boson Hamiltonian Hb is pro- 
posed as a model Hamiltonian to study properties of liq- 
uid ^He 1^ , granular superconducting arrays , and fiux 
lines in superconductors The ground state phase di- 
agram has been studied using mean-field theories ||, ||] 
and numerical methods |^ , and shows a Mott insulating 
phase, a superfluid phase, and a phase with having both 
the orders simultaneously {supersolid phase). The corre- 
spondence of these states to those in the spin model is as 
follows: Mott insulating and Neel states, superfluid and 
spin-fiopping states, and supersolid and "intermediate" 
spin states j^, respectively. 

The main purpose of this paper is to show a clear ev- 
idence for the first-order phase transition of the 2D and 
3D spin-1/2 XXZ models in the presence of a magnetic 
field by using recently developed numerical method and 
present the ground state phase diagrams. The presence of 
a phase-separated phase and the absence of the supersolid 
phase in the closely related system of hard core bosons 
with nearest-neighbor repulsions are also discussed. 

The magnetization process of the 2D and 3D spin-1/2 
XXZ models defined by Eq. (^ is studied numerically 
on square (number of spins Ns = L x L) and cubic 
{Ng = L X L X L) lattices using quantum Monte Carlo 
(QMC) technique based on stochastic series expansions 
(SSE) pH ]. Very recently an important technical im- 
provement have been achieved by Sandvik jl^ . He found 
a new algorithm of cluster-type updates (operator-loop 
updates) within the SSE QMC scheme which reduces au- 
tocorrelation time drastically compared to simulations 
using only local updates. While this method is very sim- 
ilar to the loop algorithm in the world-line QMC method 
proposed by Evertz et al. Q , one major advantage of the 
SSE method with operator-loop updates is that there is 
no difficulty in simulating systems with anisotropic cou- 
plings in external magnetic fields owing to not having 
to have "freezing" configurations and "global" weights 
which make the loop algorithm in the world-line QMC 
method highly inefficient . This reduced autocorrela- 
tion time enables us to go down to very low temperatures 
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in very high magnetic fields and therefore the method 
is suitable for our purpose. Another advantage of this 
study using the SSE scheme over other earlier numerical 
studies Q |l^, 1^ is that the simulations are performed 
directly in the ground canonical ensemble, i.e., magne- 
tization per site ruz — Mz/Ns is calculated for a given 
magnetic field h. In this paper temperatures T are set 
to be J/k-QT = 2L which is low enough to study the 
ground state properties on finite lattices and peri- 
odic boundary condition is used. The exchange coupling 
J will be taken as the energy unit. Since there exists 
long-range Neel-ordered state only for A > 1.0, from 
which the first-order spin-flop transition can take place 
by applying a flnite external field, our main focus in this 
paper will be put on this anisotropic regime. 
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FIG. 1: (a) Magnetization curve of the 2D spin-1/2 XXZ 
model with A = 1.5. Inset: enlarged scale is used. Error 
bars are smaller than symbols, (b) Correlation time t(1) (for 
definition see in the text) of (circles) and ^(Tr, tt) (squares) 
as a function of bin length I. The parameters are Ns = 12^^, 
T = 1/24, A = 1.5, and h = 2.0. Solid (open) marks are 
data for operator-loop updates with closed loops constructed 
as many as 25 (100) times per MC step, (c) The ground 
state phase diagram of the 2D spin-1/2 XXZ model with 
anisotropic constant A in the presence of the magnetic field 
h. There exist three phases: (i) Neel ordered phase, (ii) spin- 
fiopping phase, and (iii) fully saturated ferromagnetic phase. 
Solid line is 2(1 + A) (see in the text). Dashed line is a guide 
to the eye. 

Let us first study the magnetization process in 2D. A 
typical example of the magnetization curves is shown in 
Fig. (a). One can see that at certain critical magnetic 
field he the magnetization changes discontinuously 
from to a finite value . For this example in the figure 
with A = 1.5 the magnetization jumps at he 1.83 to 
'--^ 0.11 |l^. Working on various system sizes (see in 
the inset of Fig. |l|(a)), we confirm that finite size effects 
are small and conclude that the jump in magnetization 
is not due to extrinsic factors of working on finite-size 
lattices. The result already gives a clear evidence for the 
first-order phase transition. 

In order to illustrate that autocorrelation times of our 



QMC measurements are short enough, the integrated au- 
tocorrelation time Tint is estimated as follows: we first 
divide a sequence of Monte Carlo data points into bins 
of length /, and for each bin length I, the average of the 
data in the 6th bin and the variance cr^(/) of the bin 
averages are calculated. The integrated autocorrelation 
time Tint is then estimated by the asymptotic value of 
t{1) = /g^(/)/2g^(/ = 1) at large / above which t{1) does 
not depend on I pTj l . The results of Tint for and spin 
structure factor 5(q) at q = (7r,7r) are shown in Fig. |l| 
(b) for A = 1.5 and h = 2.0. From this figure Tint's are 
estimated to be no longer than 4-5 MC steps. Since we 
use bin length I — 1000 or more to estimate the statisti- 
cal errors, we can be sure that our error bar is accurate. 
It should be noted here that Tint can be controlled by 
changing how many times closed loops are constructed 
in the operator- loop updates. In this paper we have to 
construct closed loops as many as about 1000 times per 
MC step for larger values of A. 

Repeating the same procedure for different values of A, 
the ground state phase diagram of the 2D spin-1/2 XXZ 
model is completed. The result is given in Fig. Q (c). 
There exist three different phases denoted by (i) rUz — 0, 
(ii) niz ^ 0, and (iii) = 0.5 in the Figure. The third 
phase (iii) corresponds to a fully saturated ferromagnetic 
phase and is separated from the second phase (ii) by the 
critical magnetic field h™'^^. This transition is trivial and 
is not our interest here. The critical magnetic field h™^^ 
is indeed easily calculated by going to the boson model 
described by Eq. (^): the critical chemical potential to 
have just one particle in the d-dimensional boson system 
is -2dt and therefore h'^''''/J is found to be d{l + A). 

To show further evidence for the first-order transi- 
tion between phase (i) and phase (ii) we calculate the 
spin structure factor 5(q) = l/iV^ J]; j e*i ('-j)Si • Sj at 
q = (tt, tt) and the spin stiffness (helicity modulus) ps as 
a function of h for a fixed A. The spin stiffness ps is cal- 
culated by, for example, ps = T{wl+Wy{+wl)) /dL'^~^ in 
2D (3D) where Wa is the winding number per linear spa- 
tial lattice size L in the a-direction [Q. ps corresponds 
to the superfluid density in the boson model and is used 
to detect the superfluid phase in the system ||, |l^, p^ . 
The results are shown in Fig. ||(a) for A — 1.5 as a func- 
tion of h. When h is small S{q) has a peak at q = (tt, tt) 
(q dependence is not shown) and ps = 0. With increasing 
h these quantities change discontinuously at he where the 
magnetization to^ jumps, and 5'(7r, tt) becomes zero while 
Ps has a finite value ||l^. Apparently these two phases 
are (i) Neel ordered phase (for h < he) and (ii) spin- 
fiopping phase (for h > he), and are separated through 
the first-order transition pO| . 

Earlier studies found the supersolid phase in the boson 
model characterized by having finite values of S'(7r, tt) and 
Ps simultaneously in a region between phase (i) and phase 
(ii) of the phase diagram H, ^ . One can indeed barely 
see small but finite values of ^(Tr, tt) for h > he where ps is 
finite (see Fig|2| (a)). In order to elucidate systematically 
the possibility of the existence of the supersolid phase 
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in 2D, we carry out the finite size scaling analyses of 
S{tt, tt) and Ps for fixed magnetizations. The results are 
presented in Fig. §(b). One can see that S{tt,tt)/Ns (ps) 
stays finite in the limit of TVs ^ oo for h < he (h > 
he) while S{7r,7T)/Ns approaches to zero at Ng — > oo for 
h > he. Doing the similar analyses for different values of 
A and to^ it is found that whenever h > he S{ti,'k)/Ns 
goes to zero in the thermodynamic limit. We therefore 
conclude that the supersolid phase does not exist in the 
2D spin-1/2 XX Z model. The results are consistent with 
very recent studies by Batrouni and Scalettar p3. 
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FIG. 2: (a) S(7r,7r) and ps for the 2D spin-1/2 XXZ model 
with A = 1.5 and Na = 20^ as a function of the magnetic 
field h. Inset: enlarged scale is used, (b) Finite size scaling of 
S{n,n)/Ns and ps for different magnetizations with A = 1.5. 
The top figure shows S{T7,Tr)/Ns for iriz = and the mid- 
dle (bottom) figure shows ps {S{Ti,n)/Ns) for — 0.113. 
Dashed lines are liner fitting curves of QMC data, extrapo- 
lating to 1/Ns ^ 0. 

We now carry out the similar calculations for the 3D 
spin-1/2 XXZ model to study the nature of the phase 
transition induced by the external magnet field. Our 
strong evidence for the first-order transition is provided 
in Fig. |^(a). In the Figure magnetization to^, spin struc- 
ture factor 5'(7r, TT, tt), and spin stiffness ps for A = 1.5 
are plotted as a function of the external magnetic field h. 
It is clearly seen that as in the case of the 2D model 
those quantities change discontinuously at the critical 
magnetic field he- Working with different values of A, 
the ground state phase diagram is constructed and the 
result is shown in Fig. ^(b). As in 2D the phase dia- 
gram consists of three different phases: (i) Neel ordered 
phase, (ii) spin-flopping phase, and (iii) fully saturated 
ferromagnetic phase ||2^. The critical magnetic field he 
in 3D is observed to be larger compared to that in 2D 
for a given A. This is simply because the increased coor- 



dination number in 3D makes the needed magnetic field 
larger to destroy the Neel state. The similar finite size 
scaling analyses which was done for the 2D model do not 
find the supersolid phase in the 3D model. 
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FIG. 3: (a) m^, ^(Tr, tt, tt), and ps for the 3D spin-1/2 XXZ 
model with A = 1.5 and Ns = 8'^ as a function of the mag- 
netic field h. Inset: enlarged scale is used, (b) The ground 
state phase diagram of the 3D spin-1/2 XXZ model with 
anisotropic constant A in the presence of the magnetic field 
h. The diagram consists of (i) Neel ordered phase, (ii) spin- 
flopping phase, and (iii) fully saturated ferromagnetic phase. 
Sofid line is 3(1 -I- A) (see in the text). Dashed line is a guide 
to the eye. 

Finally we summarize our results in Fig. ^ by showing 
the ground state phase diagrams of the 2D and 3D spin- 
1/2 XXZ models on the parameter space of anisotropic 
constant A and magnetization . The phase diagrams 
show (i) the Neel ordered state at A > 1.0 and rriz = 0.0, 
(ii) the spin-flopping state with p^ ^ 0, (iii) the fully 
saturated ferromagnetic state at = 0.5, and (iv) a 
phase-separated state (denoted by PS in the Figures). 
The phase-separated region exists because as was seen in 
Figs, ^a) and |^(a) there is a region in magnetization tUz 
to which we can not reach in a thermodynamically stable 
way no matter how fine the magnetic field h is tuned. In 
other words, if one could have a states with rriz which is 
in this magnetization region, the state would be phase- 
separated between the Neel state with — and the 
spin- flopping state with rriz = rn'i,. Some earlier stud- 
ies predicted two phase coexisting regions working with 
the canonical ensemble (i.e., fixed m^) and their results 
were interpreted as the supersolid phase ^ . However 
our calculations conclude that these phases are thermo- 
dynamically unstable and are phase-separated. Our con- 
clusions are consistent with recent studies by Batrouni 
and Scalettar for the 2D boson model [ p^ . 

In conclusion, we have studied numerically the magne- 
tization process of the spin-1/2 anisotropic XXZ model 
in two and three spatial dimensions using QMC method 
based on SSE and shown clear evidences of the first-order 
phase transition in the presence of an external magnetic 
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FIG. 4: The ground state phase diagrams of the (a) 2D and 
(b) 3D spin-1/2 XXZ models on the parameter space of 
anisotropic constant A and magnetization m^. PS stands for 
phase separation. The Neel (fuUy saturated ferromagnetic) 
sate exists in the region of A > 1.0 and = 0.0 (m^ = 0.5). 
The spin-flopping phase is in the region denoted by ps 7^ 0. 
Dashed line is a guide to the eye. For the boson model de- 
fined by Eq. (^) and A correspond to 0.5 — n and V/2t, 
respectively. 



field. Based on the calculated ground state phase dia- 
grams, the existence of a phase-separated region and the 
absence of supersolid phases are pointed out in the re- 
lated systems of hard core bosons with nearest-neighbor 
repulsions. 



It would be of great interest to study effects of ran- 
dom anisotropic constant A on the phase diagrams 
since one can elucidate quantum effects on recently pro- 
posed impurity-induced quantum-critical-point-like be- 
havior near first-order phase transitions in the Ising 
models Q. Another interesting issue to address is to 
see effects of long-ranged Coulomb repulsions between 
bosons on the phase-separated region in the phase di- 
agram. There is a general belief that introducing the 
Coulomb interactions replaces the phase-separated state 
by thermodynamically stable states such as a droplet-like 
state and perhaps a stripe state Q . The model studied 
here provide an ideal system to examine the possibilities 
of those exotic states using unbiased numerical methods. 



The author thanks E. Dagotto and A. Dorneich for 
valuable discussions. 
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